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The Seventh Annual 
International Conference on Technology 
in Collegiate Mathematics 

November 1 7-20, 1 994 


Hosted by Valencia Community College, Sponsored by Adison-Wesley Publishing Company 
Founded by Franklin Demana and Bernd Waits, Conference Chair Judith Jones 

Walt Disney World Dolphin Hotel, Orlando Florida 

Mail or Fax to: Addison-Wesley Publishing Company, Attn: Beth Sheehan, 1 Jacob Way, 

Reading, MA 01867 Fax # 617-944-8964 


We have a FAX now!! We have a FAX now!! We have a FAX now!! 
Our FAX number is our phone number: 

+43 ( 0)2275 8207 

We have a FAX now!! We have a FAX now!! We have a FAX now!! 


At this occasion a FAX number correction: 

Herbert Appel's FAX number is ++(0)9264 8700 (Germany) 
Sorry Herbert for the mistake 
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Liebes DUG - Mitglied, 


Dear DUG Member, 


ich freue mich wirklich, Sie alle wieder in der DUG 
begriiBen zu diirfen. Bei dieser Gelegenheit mochte 
ich mich fur die vielen lieben Briefe und fur Ihre 
Anerkennung fur unsere Arbeit recht herzlich 
bedanken. Eine kleine Auswahl finden Sie im User 
Forum. Verstehen Sie das bitte nicht als Eigenlob, 
denn Ihre Zustimmung gilt vor allem den recht 
zahlreichen DERIVIANERN, die immer wieder 
Anregungen und Beitrage liefern. Ein Blick auf die 
Themenvorschau laBt ein interessantes und 
reichhaltiges DNL-Jahr erwarten. Leider laBt es 
meine Zeit (noch) nicht zu, den Umfang des DNL 
zu erweitern. (Ich habe auch noch einen 
"Nebenberuf’!) 

Ich mochte besonders auf die zuletzt versprochenen 
Algebra- und Zahlentheoriebeitrage von Dr. 
Wiesenbauer hinweisen. Dr. Wiesenbauer hat er- 
freulicherweise sein Versprechen in diesem DNL 
bereits eingelost. Ich kann keinen der vor uns lie- 
genden Beitrage hervorheben, jeder ist fur sich ori- 
ginell und einmalig. Die Auszuge aus dem BBS 
werden fortgesetzt. Ein paar Seiten Computerlisting 
(ca 60) liegen noch zur Auswertung bereit. Ich kann 
Ihnen einige interessante elektronische Diskus- 
sionen ankiindigen. 

Vor einer Woche erhielt ich eine wunderschone 
DERIVE-Animation aus Norwegen. Halvor Devoid 
erzeugte eine Show "Kepler's World System in 
3-D" mit DERIVE und GRASP. Herzlichen Dank, 
lieber Halvor. Die Teilnehmer von Krems 93 wer- 
den sich geme an Deine beeindruckende Demon- 
stration erinnem. Da Du mir die Weitergabe freige- 
stellt hast, werde ich Deine Show in die Jahresdis- 
kette 1994 aufnehmen. Ich bin sicher, damit in 
Deinem Sinn zu handeln. 

Zum SchluB noch drei administrative Hinweise: 

Wir sind endlich mit einem FAX erreichbar 
(Siehe Informations). 

Sie erhalten mit dem nachsten DNL eine Quit- 
tung iiber Ihren Mitgliedsbeitrag 1994. 

Falls Sie den Jahresbeitrag noch nicht geleistet 
haben, dann beniitzen Sie bitte das beiliegende 
Erinnerungsschreiben. Danke im voraus. 

Ich freue mich auf ein DUG Jahr 4 mit Ihnen 



Em really glad to meet you all again in the DUG. I 
would like to take the opportunity to thank for all 
your nice letters and for your appreciation. You will 
find a little selection in the User Forum. Please 
don't see this as self-praise, because your approval 
especially goes to the numerous DERIVIANS, who 
submit papers and give suggestions. One glance at 
the preview shows that we can expect an interesting 
and contentful DNL-year. I'm sorry to say that I 
don't have time (not yet) to extend the content of 
the DNL (I have a "sideline", too!) 

I want to draw your attention to Dr Wiesenbauer 's 
Algebra and Number Theory contributions promi- 
sed in the last DNL. I'm very glad that Dr Wiesen- 
bauer has fulfilled his promise. It is really impossi- 
ble for me to particularly point out one or the other 
contribution in the next issues, each one of them is 
original and unique. The extracts from the BBS will 
be continued. Some pages of computer paper ( « 60) 
are waiting to be evaluated, whether they could be 
of interest for you. I can promise some useful 
electronic discussions 

One week ago I received a wonderful DERIVE 
animation from Norway. Halvor Devoid produced a 
show "Kepler's World System in 3-D". Many 
thanks, dear Halvor. Participants of the Krems 
Conference 1993 will certainly remember your im- 
pressive demonstration. As you have given your 
permission to distribute the file I feel free to 
incorporate your show into the 1994 diskette at the 
end of the year. I'm sure to act in your sense. 
Thanks again. 

Three administrative notes at last: 

You can reach us by FAX now . If it does not 
work the first time, please give it another try. 
(See Informationpage) 

Togetherwith the next DNL you will receive a 
receipt about your membership fees for 1994. 

If you haven't settled your dues for 1994, then 
please use the reminder. Thanks in advance. 

I am looking forward to the 4th DUG year with you 

Sincerely yours 



Ihr 





P 2 


EDITORIAL 


D-N-L#13 


The DERIVE-NEWSLETTER is the 
Bulletin of the DERIVE User Group. It is 
published at least four times a year with a 
contents of 30 pages minimum. The goals 
of the D-N-L are to enable the exchange of 
experiences made with DERIVE as well as 
to create a group to discuss the possi- 
bilities of new methodical and didactical 
manners in teaching mathematics. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/Fax: +43-(0)2275/8207 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to 
reinforce the international touch of the 
D-N-L. It must be said, though, that non- 
English articles will be warmly welcomed 
nonetheless. Your contributions will be 
edited but not assessed. By submitting 
articles the author gives his permission for 
reprinting it in D-N-L. The more contri- 
butions you will send, the more lively and 
richer in contents the DERIVE Newsletter 
will be. 


Preview: (Contributions for the next issues): 

Fluid flow in DERIVE, Reuthera.o., BRA 

Applications in Electrical Engineering, Scheuermann, GER 

Stability od systems of ODEs, Kozubik, SLK 

Los desplazamientos en las finciones elementales, Ramos, ESP 

Algebraic Operations on Polynomials in DERIVE, Roanes, ESP 

Prime Iterating Number Generators, Wild, UK 

Graphic Integration, Probability Theory, Linear Programming, Bohm, AUS 
DERIVE in Austrian Schools, some examples, Lechner & Eisler, AUS 
Tilgung fremderregter Schwingungen, Klingen, GER, 

Der Fermat-Punkt im Dreieck, Geyer, GER 
Continued Fractions in DERIVE, Cordoba a.o., ESP 
Life Game, Turtle-Commands in DERIVE, Lechner, AUS 
Dreieck. MTH, Wadsack, AUS, and others 

At this position in the last DNL you could find my Christmas-tree. In the letter of the Editor I 
expressed my hopes that I wouldn't have to draw a DERIVE Easter Bunny. But look, Mr 
Klingen from Bonn sent a letter containing a wonderful rabbit: 


....Because you have wished an Easter Bunny 
in your last issue, here is one. It is built from 
28 points with 5 parametric splines 
consisting each of 6 or 7 points. 



Impressum: 

Medieninhaber : DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber : Mag . Josef Bohm 
Herstellung: Selbstverlag 
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(Menemui is the abbreviation for Menemui Matematik published by the Malaysian Mathematical Society.) 

Ob diese Aufgabe "zu FuB" gelost werden kann? Mit DERIVE geht es in einer Zeile (allerdings wird 
MISC.MTH gebraucht) und hat bei mir mit DERIVE XM etwa 6,5 Minuten gedauert. Tatsachlich hatte ich die 
Losung zunachst durch Expandieren, Speichem auf eine Datei und " Anschauen" gefunden. 


1: "MISC.MTH loaded with LOAD UTILITY" 

2 : MAX (VECTOR (POLY_COEFF ( ( 1+2 *x+3*x A 2+4 *x A 3 ) A 20,x,i),i,60)) 

3: 88423110877597693745 

4. " 387,2 sec on a 486/33MHz with DERIVE XM 

Is there anybody who knows how to solve this problem by reasoning? My first attempt was to EXPAND, to 
SAVE the expression and look for the largest coefficient. 


Derive 6 needs only 2 seconds (without loading the utility file). 


Paul Drijvers, Utrecht, Netherlands 

Can you draw this picture in DERIVE? (Do not look at the back! 


DNL: I did. And thanks, Paul for your nice tree: 


f 1 15 ) 

#1: F(t) := 2-t + + 1 (10 t - FLOOP(lOt)) 

K 2 10 J 

#2: G(t) := 10 - 10- 1 


#3: 

#4: 


#5: 


1 

H(t) - + t 

2 

k(t) := 5ASIN(SIN(nt)) 
F(t) G(t) 

- F(t) G(t) 

. H(t) k(t) . 


[TJ^HPlot tt 5 with t from 0 to 1 



Dr. Wolfram Koepf, Berlin, Germany 

Josef Lechner wants a fifth argument of the ITERATES function to construct, e.g., the expression 

Do A {xJ),f 2 M (x 0 )),f 3 (f 2 (f, (x 0 )))....] 

(Bulletin of the DERIVE User Group 12, page 4). I agree absolutely with this wish! However, I like to show, 
how one is able to solve the given problem with the current versions of DERIVE: recursively. 

The DERIVE function: 


GENERATE (list, k) := 

IF (k=l, [ELEMENT (list, 1) ] , APPEND (GENERATE (list, k-1) , 

[LIM (ELEMENT (list,k) ,x, ELEMENT (GENERATE (list, k-1) ,k-l) ) ] ) ) 

does the Job, required. If we declare F1,F2,F3 and F4 as functions of one variable, too, and define 
list : = [FI (x) , F2 (x) , F3 (x) , F4 (x) ] , then simplification of the expression 
GENERATE (list, DIMENSION (list) ) leads to the result 

InputNode := Word 

[fl(x) f2(x) f3(x) := f400 :=] 

GENEFATE([ f 1 (x) , f2(x), f4(x>], 4) 

[fl(x), f2(fl 00), f3Cf2CflCx»), f4(f3Cf2CflCx»)>] 



D 


-N-L#13 


DERIVE 


USER 


FORUM 


P 5 


If we declare, however F5, F10 as functions of one variable, too, and redefine 

list : = [FI (x) , F2 (x) , F3 (x) , F4 (x) , F5 (x) , F6 (x) , F7 (x) , F8 (x) , F9 (x) , F10 (x) ], then 
simplification of GENERATE (list, DIMENSION (list) ) leads to memory overflow: In each step the 
recursion invokes the previous value twice, which generates hyperexponential complexity, and hyperexponential 
memory requirements. This problem can be resolved using an auxiliary function. We redefine GENERATE by 
the statements 


GENERATE_AUX (list, k, aux) := 

APPEND (aux, [LIM (ELEMENT (list,k) , x, ELEMENT (aux,k-l) ) ] ) 

GENERATE (list, k) := 

IF (k=l , [ELEMENT (list, 1) ] , GENERATE_AUX ( list , k, GENERATE (list, k-1) ) ) 


and simplification of the expression GENERATE (list, DIMENSION (list) ) leads to 


[f5(x) :=, f6(x) :=, f7(x) :=, fS(x) :=, f9(x) :=, flO(x) :=] 

GENERATE([fl(x) , f2(x), f3(x), f4(x), f5(x), f6(x) , f7(x), fS(x), f9(x) , flO(x)], 10) 

[fl(x), f2(fl (x)) s f3 (f2(fl (x))) , f4 (f 3 (f 2 (f 1 (x) ) ) ) , f5(f4(f3(f2(fl(x))))) , f6(f5(f4(f3(f2(fl(x)))))) , 
f7 (f6(f5 (f4(f3 (f2(fl (x))))))) , fS(f7 (f6(f5 (f4(f3 (f2(fl (x)))))))) , 
f9CfS(f7Cf6(f5Cf4Cf3Cf2Cfl(x))))))))) ! fl0Cf9CfSCf7(f6Cf5Cf4Cf3Cf2(flCx))))))))))] 


I checked, but in case of 10 nested functions the auxiliary version is not necessary now. See some 
examples: 

This is the example given in DNL#1 2 



2 

x 1 


GENERATE 

SIN(x) , 1 + x, , 

. 4 


2 x . 



2 



SIN(x) 



SIN(x), SIN(x) + 1, + SIN(x) + , 

2 2 2 

SIN(x) + 2*SIN(x) + 1 


list 



1 - 2 + x, 


1 


x 


x 3 

GENERATE(1 ist , DIHENSION(list)) 


2 

+ 3*x, 1 + 4*x, 1 + 5*x, 1 + 6+x, x - 1 


2 

2 2 1 1 2 + (x - 1) 4 20 

x, x i 1 — 2 + x i ) j j 5 — i 26 — i 157 — 

2 2 2 2 2 
1 - 2*x 3 ■ (1 - 2.x ) 2.x - 1 2.x - 1 2*x - 1 

2 

120 240. (314.x - 217) 

, 24643 

2 2 2 
2.x - 1 (2.x - 1) 


A.C. Robin, Colchester, England 

I note that my difficulty with defining the derivative of x A 3, using Derive has provoked some interest. Dr 
Schumm's solution using 'LIMIT ' in DNL#12 is a more satisfactory solution than the use of DCUBE(u,x) in 
DNL#11, as whilst this gave DCUBE(x+z,x) correctly it did not give DCUBE(x A 2+3,x) as 3(x A 2+3) A 2 as, of 
course I wished. 
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I have another similar problem involving integrals, and the integrals I have been using often involve 
INT(e A (-l) * x A n,x,0,inf). This integral should simplify to n!, Derive does not do this simplification, but it does 
verify it, for actual values of n. What I wish to know is whether there is any way of teaching Derive how to do 
integrals with which it is not familiar (thinking of this one in particular), so that when it encounters it in a more 
complicated expression, it can be simplified correctly. 


PS.: Has anyone thought of having a DERIVE Bulletin Board in the UK? 

J.M.M.Cardia Lopes, Porto, Portugal 

I have tried to run DERIVE 2.03 at a HP 9000 750/33 workstation by means of SOFT-PC (a MS-DOS emulator 
for UNIX). DERIVE works but the first screen stays over the next screens. The same with the graphics 
generated by PLOT. It looks bad and is very deficient. 

Do you know how to solve this bug? Do you know anyone who has the same problem? 


K. Winkelhausen, Arch, Switzerland 

Ich ware froh, wenn ich einen Ausdruck samtlicher DERIVE- Anwender der Schweiz erhalten konnte, um einen 
regen Kontaktaustausch innerhalb des Landes vomehmen zu konnen. 

DNL: Lieber Herr Winkelhausen! Ich mochte grundsatzlich keine Adressen weiterleiten. Ich nehme 
aber an, daft Sie mit der Angabe Ihrer Anschrift einverstanden sind, so konnen sich alle Schweizer 
DERIVE-Anwender mit Ihnen in Verbindung setzen: Solothurnstr. 15, CH-3296, Arch. 

K. Herdt, Osnabriick, Germany 

Fur die Einfuhrung und Erlauterung der Darstellung von Funktionen in Polarkoordinaten fmde ich es aus 
didaktischen Griinden sehr hilfreich, daB DERIVE auch hierbei eine (undokumentierte?!) Parameterdarstellung 
erlaubt. Bei der Einstellung des Darstellungsmodus ’’Polar" wird namlich der erste Term des geordneten Paares 
[a(t),b(t)] als Parameterdarstellung des Leitstrahls/Radiusvektors r interpretiert, der zweite als Parameter-dar- 
stellung des Polarwinkels cp. So wird beispielsweise durch [t, 1] eine unter dem Winkel 1 (rad!) vom Ursprung 
ausgehende Strecke dargestellt, deren Lange durch den Laufbereich des Parameters t gegeben ist. 

Demgegeniiber wird [l,t] als Kreisbogenstuck vom Radius 1 interpretiert und [t,t] naturlich als Spirale. 

Soil nun beispielsweise fur die in Polarkoordinaten gegebene Funktion r = sin(2cp) eine Unterteilung des 
Defmitionsbereichs [0,27i] in vier Teilintervalle erfolgen, so laBt sich dies erreichen durch: 

VECTOR ( [SIN (2 (k*7i/2+t) ) ,k*:i/2+t] ,k, 0, 3) 

anschlieBend Simplify und Plot, wobei der Parameter tjeweils von 0 bis 7i/2 laufen muB (einmalige Eingabe mit 
Ctrl-Enter abschlieBen). Sofern die Farbeinstellung fur die Darstellung der Graphen auf Auto steht, werden nun 
die vier Teilabschnitte des "Vierblatts" in unterschiedlichen Farben dargestellt und die Zuordnung dieser Teile 
zu den Abschnitten des Defmitionsbereichs laBt sich besser erkennen. 

DERIVE allows an - undocumented? - 
parameter form for functions in polar form. This 
is very helpful for introduction and explication of 
functions in polar form. Having set the 
coordinate system to polar the expression 
[a(t), b(t)] is interpreted as parameter form 
[r(t), cp(t)]. So [t, 1] is a segment of legth t which 
forms an angle of 1 rad with the x-axis. [1 , t] is a 
part of the circumference of a circle with radius 
1 and [t, t] is a spiral. 

Expression #1 gives a demonstration of a 
quadrifolium r = sin(2cp) (in polar form) which 
shows the four parts of the leave in different 
colours (Option > Display > Color > 

Automatically 


change color of new plots) and makes clear the 
assignment of these parts to the parts of the 
domain 


0 <t<nl2 


/I : VECTOR^ SIN| 2 

0 < t < jt/1 2 
* 2 \ VECTOr| 

0< t< nl2 


(H 2 '("T’ ' '))’ ~T~ * *]’ "• ’) 


k ♦ 7T 

1, + t 

6 


, k, 0, 12^ 


#3 : VECTOR 


k » 7T k ♦ 7T 

+ t, + t 

2 2 


, k, 0, 7^ 
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Weniger erfreulich sind die von Version zu Version vorgenommenen (undokumentierten !) kleineren 
Anderungen bei bereits vorhandenen Funktionen und Operatoren von DERIVE. Bis zur Version 2.58 wurde 
streng unterschieden zwischen den drei Objekten 

Vektor (= geordnete Liste von Elementen), 
einzeilige Matrix und 

einspaltige Matrix (vergl. auch Handbuch, Vers. 2, S.135 ff). 

Verstandlicherweise lieB sich daher ein Vektor (im Unterschied zu den beiden anderen Objekten) nicht 
transponieren. Ab Version 2.59 laBt sich nun plotzlich auch ein Vektor transponieren und es entsteht eine 
einspaltige Matrix, die konsequenterweise bei emeutem Transponieren in eine einzeilige Matrix iibergeht. Fur 
jeden vom Nullvektor verschiedenen Vektor v ergibt sich also die sicherlich nicht sonderlich befriedigende 
Ungleichung: . 

Da die Autoren von DERIVE auch mit einem bestimmten padagogischen Anspruch angetreten sind, ist dies 
eigentlich bedauerlich. 

DNL : Answers to inquiries of Mr Herdt and Mr Lopes: 

In one of the later messages in the BBS I found the similar problem: 

Message 3374 From SOFT WAREHOUSE to KEITH WILLIAMS about TRANPOSE OPERATOR BUG 

Thank you for your fax letter dated December 24, 1993 concerning a problem with DERIVE’s transpose 
operator. The problem arose because DERIVE version 2.59 improperly transposes ragged arrays. Ragged arrays 
are vectors of vectors not all having the same number of elements. The problem did not arise in earlier versions 
of DERIVE because those versions did not even attempt to transpose ragged arrays. 

The problem has been resolved and will be included in the next version of DERIVE. Please send me your 
mailing address and I will send you a complementary update as soon as it is available. Aloha, A1 Rich, 

Soft Warehouse, Inc. 


I sent a FAX to SWHH in the evening and thanks modern times same night at 03:22 am I received an 
answer from Al Rich: 

In response to requests by users, beginning with version 2.59 DERIVE transposes a vector to an n by 1 column 
matrix. The transpose of an n by 1 column matrix is a n by 1 row matrix. Since a 1 by n row matrix is not 
equivalent to an n element vector they cannot be subtracted and their difference is not the zero vector. 

Perhaps DERIVE should transform the vector into a row matrix so the subtraction can occur. However, this is 
not valid if the vector has symbolic elements that couls themselves be replaced by vectors. 


The screen problem that occurs when running DERIVE under Soft-PC must be a bug in Soft-PC, not DERIVE. 
Apparently Soft-PC does not correctly emulate the screen services provided by the IBM-PC BIOS. I recommend 
Mr. Lopes contact the authors of Soft-PC to see if they can resolve the problem. Sicerely, Al Rich 


p 
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Message 2392: From SOFT WAREHOUSE to PUBLIC about L'HOPITAL'S RULE 

The Derive LIM function is quite powerful, but it is educational to implement your own function for computing 
limits using L’Hopital's rule particularly if you want the rule to be automatically repeated if necessary. 
Consequently, file LHOPITAL.MTH is an interactive lesson that helps teach recursive Derive programming. 

To take the lesson, begin with Transfer Load Derive LHOPITAL, then follow the directions on the screen. Note 
that the use of LIM within the definition in the file is merely to SUBSTITUTE the limit abscissa for the limit 
variable in the numerator and in the denominator, so it is not cheating by using the internal implementation of 
L'Hopital's rule. 


File LHOPITAL. DFW: Recursive L'Hopital's rule for 0/0 

L(n, d, x , a) := 

If LIH(n , x, a) = 0 A Lim(d, x, a) = 0 
#1: L(0(n , x), d(d, x) , x, a) 

LUlCn, x, a)/Lm(d, x, a) 


Simplify the following two examples, & note how DERIVE returns a conditional expression when it 
cannot determine that a predicate is always true or always false. 

*2: L(5IN(x), 5 INH(x), x, 0) 

#3: L(5IN(a + x) , SINH(b*x), x, 0) 

Then, use F2 and F3 to redefine L to have a 4th 'unknown' clause identical to its else clause, 
and resimplify the second example. 

First Example: 

#4: L(SIN(x), SINHCx), x, 0) = 1 

Second Example: 

#5: L(SIN(a.x), SINH(b.x), x, 0) 


# 6 : 


IF 


a = 0 A b = 0, L 


d d 

— (a.COS(a.x)), — 
dx dx 


b • x - b • x 

b *c b*c 

+ 


» x, 0 


1 im 
x-»0 


lim a*COS(a*x) 
x-»0 


Highlight the third expression and simplify #6 once more: 




r 


' b • x - b • x " 






d d 

b*c b*c 


a 

#7: 

IF 

a = 0 A b = 0, L 

— (a.COS(a.x)), — 

+ 

, x, 0 

i 




.dx dx 

L 2 2 J 


b J 


1 im a*COS(a»x) = a 
#8: x-0 


#9: 


( b • x 

b*c 

1 im 

x->0 L 2 


+ 


- b*x ) 

b*c 


2 


= b 


The recommended redefinition of L(n,d,x,a): 

L(n, d, x, a) := 

If LUn(n, x, a) = 0 a LIH(d, x, a) = 0 
*10: L(a(n, x), 3(d, x), x, a) 

Lm{n, x, a)/LU*l(d, x, a) 

LU1(n. x, a)/LH>l(d, x, a) 


*11: L(SIN(a.x), SINH(b.x), x, 0) 



I recommend using DERIVE's stepwise simplification feature. Josef 
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Message 2397: From SOFT WAREHOUSE to PUBLIC about MOCK CHAOS 

When plotted using the default range -pi < t < pi, the following elementary expression closely mimics the usual 
phase plot of the chaotic Lorenz strange attractor!: 


[ATAN (t) + t SIN (t A 3) , ATAN (t) + COS (t A 3)] 



Message 2399: From SOFT WAREHOUSE to PUBLIC about A DERIVE DATA-BASE PROGRAM 

File DATABASE.MTH implements a simple data base program and gives examples of its use. 

"File DataBase . MTH : Example of data-base table look-up:" 

"1st establish some unassigned multi-character vars even in character mode:" 

age : - 

height : = 

eye_color : 

b r own : - 

blue : - 

"Establish some data records:" 
darth_vader : = [ [age, 84 ] , [weight, 240 ] ] 

veronica : = [ [age, 30 ] , [weight, 130 ] , [eye_color, brown] ] 

"GET_AUX decrements k until it is 0 or the key is found:" 

"Note how we must use 3-valued logic (true, false, unknown) because Derive" 

"usually can't determine if two distinct nonnumeric keys are identically" 

"equal as math vars or never equal as math vars." 

"For example, the equation 'age=weight' could be identically true, or never true. 
GET_AUX (rec, key, k) :=IF (k=0, ?, IF ( key=ELEMENT (rec, k, 1) , ELEMENT (rec, k, 2) , G 

GET_AUX (rec, key, k-1) , GET_AUX ( rec, key , k-1 ) ) ) 

"Use Manage Substitute to replace d with the appropriate subexpression" 

" so that the subsequent two test cases simplify to 'brown' and '?':" 

GET (rec, key) : =GET__AUX (rec, key, d) 

GET (veronica, eye_color) 

GET (darth_vader , eye_color) 

d::3 

I tried d = 3: 


GET(veronica, eye_color) = brown 
GET(darth_vader , eye_color) = ? 
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Message 2538: From WILLIAMSON to SOFT WAREHOUSE about #2399 / DAT AB ASE.MTH 
COMMENTS 

The instructions with DATABASE. MTH could be clearer. It should be pointed out that the k parameter in 
GET AUX represents the number of fields in a record that can contain an undetermined number of fields. 
Dennis 


Message 2583: From SOFT WAREHOUSE to PUBLIC about MATRIX KRONECKER 

File KRONPROD.MTH contains a function that returns the Kronecker Product of two matrices. 


NCOLS (a) : =DIM(a;l) 

KRON_PROD_AUX (a, b, m, n, p, q) -VECTOR (VECTOR (a; ( 1+FL00R ( i , p) ) 4 . (1 + 

FLOOR ( j, q) )bi ( 1+MOD ( i , p) ) i (l+MOD(j,q) ) , j , 0 , n • q-1 ) , i , 0 , m • p-1 ) 
KRON_PROD (a,b) - KRON_PROD_AUX (a, b, DIM (a) , NCOLS (a) , DIM(b) , NCOLS (b) ) 

See one example: 


ml := 


2 3 4 0 

5 6 7 8 


KRON_PROD(ml, m2) = 


1 

2 

3 


4 

5 

6 


7 

8 

9 



2 

4 

6 

3 

6 

9 

4 

8 

12 

0 

0 

0 

8 

10 

12 

12 

15 

18 

16 

20 

24 

0 

0 

0 

14 

16 

18 

21 

24 

27 

28 

32 

36 

0 

0 

0 


5 

10 

15 

6 

12 

18 

7 

14 

21 

8 

16 

24 

20 

25 

30 

24 

30 

36 

28 

35 

42 

32 

40 

48 

35 

40 

45 

42 

48 

54 

49 

56 

63 

56 

64 

72 


A short explanation for those (few) of you who - like me - don't know the definition of the Kronecker 
product of two matrices: all elements of m2 are multiplied by the elements of ml and these products 
replace the elements of ml. The black rectangle is the product of 2 and m2 and the red one is the 
product of 7 and m2. (Definition found in the fine book of Karsten Schmidt & Gotz Trenkler, Einfuhrung 
in die Moderne Matrix Algebra, Springer.) 


Message 2607: From SOFT WAREHOUSE to PUBLIC about SORTING 

File SORT.MTH is an interactive lesson that helps teach recursive Derive programming. The goal is a function 
that sorts a vector of rational numbers into nondescending order. To take the lesson, begin with Transfer Load 
Derive SORT, then follow the directions on the screen. Your rewards will be a greater problem-solving ability 
and a useful function to add to your personal utility library. 

I renamed SORT.MTH to SORT_.MTH and the function SORT to SORT_, because SORT has been 
implemented into DERIVE since long. Josef 


"File SORT_.MTH" 

"Selection-sort elements of numeric vector v into nondescending order." 
"Find minimum element of v from element k, with element m min so far:" 
FIND_MIN (v, k, m) : =IF (k>DIM (v) , m, IF (ELEMENT (v, k) <ELEMENT (v, m) , 

FIND_MIN (v, k+1, k) , FIND_MIN (v, k+l,m) ) ) 


D 
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"Sort elements of v into nondescending order from element i on:" 
SORT_AUX ( v, i ) : =IF ( i>=DIM (v),v, ~ 

SORT_AUX (SWAP_ELEMENTS (v, i, FIND_MIN (v, i, i) ) , i+1) ) 

SORT_ ( v) := 

"Define SORT_(v), then test & debug it on the following examples: 
SORT_ ( [ ] ) 

SORT_ ( [42] ) 

SORT_ ( [5, -2, 1/2, 3] ) 


Solution and check using built-in SORT-function: 

SORT_(v) := SORT_AUX(v , 1) 


SORT_( [ ] ) = [] 


S0RT_([42]) = [42] 



1 


1 

SORT_ 

5, -2, , 3 

— 

-2, , 3, 5 


2 


2 


SORT_([77, - J2, 0, -100.25, 100- e]) 


401 

, - J2, 0, 77, 100 - e 

4 


SORTQtz, - J 2 , 0, -100.25, 100. e]) 


401 

, - J2, 0, 77, 100 * 6 

4 


Message 2699: From JERRY GLYNN to PUBLIC about OFFICIAL BUG REPORT 

FLOOR, MOD and IF functions do not work correctly in Derive with Limits. This is a known problem to A1 
and Dave who write the program but they do not have a fix for it in hand. Hopefully in time they will work it 
out. For now do not trust it. 

Message 2822: From HARALD LANG to PUBLIC about DERIVE AND EDUCATION 

There seems to be a discussion going on between Gary Ingram and Jerry Glynn, of which only parts are readable 
for the "public”, so maybe I haven’t got the issue properly. I don’t quite see why DERIVE should be viewed as a 
pedagogical tool only. I teach mathematics at the Royal Institute of Technology in Stockholm, but I use 
DERIVE there only to prepare examples and exams, but I also do research in economics at Stockholm 
University, and there I use DERIVE also for my research. To me, DERIVE is a MATHEMATICAL 
ASSISTANT, i.e., I don’t use it for anything I couldn’t in principle do on paper myself, but if I instruct my 
assistant, it does it much faster and with fewer (hopefully, no) errors. It is like a pocket calculator: it is useful 
both in education and research. To me, a useful upgrade should contain better and more efficient algorithms, 
and I agree with Gary Ingram that the latest upgrade was a bit thin in that respect. 

For instance, one thing that has irritated me is that if I solve a problem and DERIVE is unable to simplify the 
answer, I don’t get any answer at all. To fix the ideas, consider solving x from the equation ax + b = c. If a, b 
and c are messy expressions, then DERIVE tries to simplify the answer x = (c-b)/a by cancelling common 
factors in the denominator and numerator. This can easily be too complicated for DERIVE (which is natural), 
and then we don’t get any answer at all! It would be much better if DERIVE automatically gave up a 
simplification if it is too complicated, and performed a less ambitious simplification (e.g., in the example, 
simplifying the numerator and denominator separately, without trying to find common factors). Maybe this is 
the reason why MAPLE — according to Gary I, I don’t have any experience of MAPLE — is capable of solving 
more "advanced" equations? — Harald Lang 


P12 


Selections from the DERIVE-BBS 


D-N-L#13 


Message 2848: From HEINZ to PUBLIC about IN DEFENSE OF DERIVE 

I rise to speak in defense and in praise of DERIVE. I belong to that silent minority of DERIVE users who are 
not professional mathematicians. We were bom genetically deficient in math talent but with a superabundance 
of appreciation for the beauty and power of mathematics and a burning desire to solve math problems, simply 
for the joy of seeing math formulas and symbols come alive in two and three dimensional graphs. DERIVE, like 
no other program, has done that for me and for that I am very grateful. I am also grateful for the many excellent 
books that were published exploring and explaining DERIVE. Each one of these books (seven of them) opened 
up new insights into the power and and enjoyment of math and DERIVE. 

I am not completely ignorant of math. I have a PH.D. in organic chemistry and some sixty publications in the 
fields of metallo-organic complexes, catalysis and coal chemistry. But, I am 82 years old, retired and don’t have 
time to learn the language necessary to run MAPLE. I also have MathCad and like it very much, but if I want to 
take a quick look at the graph of a polynomial or do some symbolic linear algebra I turn to DERIVE. Let those 
who need a more sophisticated program use MAPLE, but don’t spoil for us, less gifted mathematically, the joy 
of using DERIVE. Heinz 


Message 2860: From JERRY GLYNN to PUBLIC about WEB FUNCTION 

On my recent trip to Seattle to the National Council of Teachers of Mathematics Meeting I had a chance to visit 
Janet Ray and Mike Pepe at Seattle Central Comm. College. They are long time users of Derive and most other 
good math software. Janet and I had discussed how to use Derive to plot some functions in Derive using the 
iterates command and this lead to discussions with Mike (this was about 1 5 months ago) and on this visit I saw 
the result of this work. It represents to me the best example of the use of Derive in a graphical manner. I strongly 
suggest that you try this out and let all of us know what you find. 


f (x) : = 

tri(x) : = [ [x,x] , [x, f (x) ] , [f(x),f(x)]] 

next (m) :=tri (element (m, 3, 1) ) 

web (a, n) : =iterates (next (m) , m, tri (a) , n) 

f (x) : =2 . 6x ( 1-x) 

web (.3,5) 

approximate the last line. Open a 2D graphing window and 
do Option State Connected and plot the result from 
approximating web (.3,5). You also should plot x and 
f(x) to see the full effect. 



I set up a list [ f ( x ) , x , web (.3,5) ] and graph it to get all 3 at once. In the function where I used 2.6 try 3.5 
or 2.1 and see the differences that emerge. As you get a feel for the situation and the time it takes your machine 
to do the calculations increase the second arguement of web. On my 486 with Derivexm I run 50 or 100 quite 
easily. 

Now you can immediately plot web(.3,5) now without approximating but activate the Option 
Approximate Before plotting in the 2D-Plot Window. 


Message 2885: From WKSARTORY to PUBLIC about COMMANDS BY FILE 

I am interested in issuing commands (such as SIMPLIFY) by putting them in a .MTH file (or some other input 
file) rather than issuing them directly through the command line interface. Is there any reference that explains 
how this can be done? «< Reply - see message #2980 »> 

Message 2971: From SOFT WAREHOUSE to PUBLIC about DERIVE MEMORY USAGE 

There has been some discussion recently on this BBS concerning memory usage by the conventional memory 
version of DERIVE. It is a complicated subject but for those that are interested perhaps I can shed some light on 
it: 

DERIVE is written in muLISP-90, which in turn is written in 8086 assembly language. The 8086 
microprocessor is a nightmare for implementing pointer-based languages like LISP. The 20 bit (one megabyte) 
address space is segmented into small 64K segments and only 4 segment registers are available. 
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Ebene Algebraische und 
T r a n s z e n d e n t e Kurven (3) 

Thomas Weth, Wurzburg, Germany 


Die Trisektrix von Maclaurin 

In der ersten Folge des DERIVE-Kurvenlexikons 
wurde die Kissoide vorgestellt. Sie diente zur 
Verdoppelung des Wurfels. Neben der Wiirfel- 
verdoppelung war die Dreiteilung eines beliebigen 
Winkels ein weiteres "klassisches" Problem. Im 
Laufe der Zeit beschaftigten sich ungezahlte 
Mathematiker und Laien damit, die 
Winkeldrittelung mit Zirkel und Lineal zu losen. 
Mit den Methoden der Algebra kann mittlerweile 
leicht gezeigt werden (vgl. Kunz), daB sich mit 
Zirkel und Lineal alleine ein beliebiger Winkel 
nicht in drei kongruente Teile zerlegen laBt. 
Trotzdem gingen und gehen immer noch 

Vorschlage von Laien an den Universitaten ein, die 
dieses Problem scheinbar elementar mit Zirkel und 
Lineal gelost haben. Unter diesen Vorschlagen 
fanden sich mitunter so komplizierte und genaue 
Naherungslosungen, daB 1983 Dudley ein 

"Handbuch" fur die gepeinigten Uberpriifer derar- 
tiger Beweise unter dem Titel "What to do, when 
the trisector comes" herausgab. 

Unter den nichtelementaren Losungen der Winkel- 
dreiteilung (also denen, welche auBer Zirkel und 
Lineal andere Hilfsmittel verwenden) linden sich 
viele, welche algebraische Kurven zu Hilfe nehmen. 
Eine der bekanntesten ist die Kurve, die von 
Maclaurin erfunden wurde. Der Schotte Colin 
Maclaurin (1698-1746) wurde mit 19 Jahren 
Professor fur Mathematik in Aberdeen, ab 1726 in 
Edinburgh. In seinem wichtigsten Werk "A treatise 
of fluxions" (1742) findet sich u.a. die nach ihm be- 
nannte Reihe sowie das gewohnlich Cauchy zuge- 
schriebene Integralkriterium fur die Konvergenz 
unendlicher Reihen. 

Auf der nachsten Seite konnen Sie die GeoGebra- 
Konstruktion sehen. AuBerdem habe ich die 
Trisektrix mit DERIVE und den Schiebereglern 
nach der Konstruktionsvorschrift hergestellt. Mit 
dem Schieberegler fur a verandert man die GroBe 
der Schlinge und der Schieberegler fur p lasst den 
Punkt P auf der Ortslinie wandem. 


Maclaurin 's Trisectrix 

in the first part of the DERIVE - lexicon of curves 
the Cissoid has been introduced. This curve was 
used for doubling a die. Beside doubling a die 
dividing an arbitrary angle into three equal parts - 
trisection - was another "classic" problem. During 
the centuries innumerable mathematicians have 
been busy to solve trisection of an angle using only 
ruler and compasses. 

By methods of algebra easily can be shown, that this 
is impossible. Nevertheless there are ideas from 
nonprof mathematicians how to solve the problem. 
There were so many very complicated ideas and 
constructions that in 1983 Dudley published a 
manual for the examiners "What do do, when 
trisector comes". 

Among the nonelementar solutions (using other 
tools beside ruler and compasses) you can find a lot 
which are using algebraic curves. One of the best 
known is the curve discovered by the Scot 
Maclaurin (1698 - 1746, professor for mathematics 
in Aberdeen with age 19, since 1726 in Edinburgh). 
In his most important opus "A treatise of fluxions" 
(1742) is to be found the socalled Maclaurin series, 
and the integral criterion for the convergence of 
infinite series (which is usually named after 
Cauchy!) 



Given are the two points S and O with SO = 2a. Two rays 
are rotating round S and O with a start position OS in the 
x-axis. The ray round S is rotating three times faster than 
the other one. The intersection P is one point of the curve 
wanted and § = ZPOS = 3 vj/ with \\r = ZPS and x-axis. 
The graphic is a Cabri-Geometre plot. 
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Konstruktion der Kurve 


Gegeben sind zwei Punkte S und O mit der Entfemung 2a. Um S und O drehen sich zwei Strahlen, beide mit 
gleichem Drehsinn und beginnend aus der Ruhelage auf der Geraden OS. Der Strahl um S dreht sich dabei 
dreimal so schnell wie der Strahl um O. Die Schnittpunkte dieser Strahlen bilden die Trisektrix von MacLaurin, 
d.h. 4>=3vi/. (Obiges Bild wurde mit Cabri- Geometre mit Hilfe der Ortslinienfunktion erstellt.) 



The left figure shows the GeoGebra-construction of the Trisectrix and the right one is the DERIVE 
construction running parallel to the analytic derivation of the equation using slider bars for a and 
parameter p. 


Die Winkeldrittelung mit Hilfe der Trisektrix 

Zunachst zeichnet man eine Trisektrix (z.B. nach 
obiger punktweiser Konstruktion). Der zu 
drittelnde Winkel wird sodann mit dem Scheitel in 
S an der x-Achse angetragen. Der freie Schenkel 
des Winkels schneidet die Trisektrix im Punkt P. 

Der Winkel ZSOP ist dann das gesuchte Drittel des 
urspriinglichen Winkels. 

Herleitung der Kurvengleichung 

Zur Herleitung der Kurvengleichung betrachtet man AOSP. Da § AuBenwinkel von AOSP ist und 4>=3v|/, gilt 
ZSPO = 2\| /. Nach dem Sinussatz gilt also in AOSP: (Sine rule!) 

r sin y sin y a a 

— = = , also r = = 

2 a sin 2 y 2 sin y cos y cosy 

cos — 

3 

Ersetzt man r durch p gemaB p sin y = r sin 3 y (Sinussatz; Sine rule) und beriicksichtigt, dass 

*2 a 3 

sin 3y = 3 sin y - 4 sin y , so erhalt man: p siny = (3siny - 4sin y), also 

cos y 

a 2 a 2 

p= (3 - 4 sin y) = (4 cos y - 1) 

cosy cosy 

Damit erhalt man die Polargleichung mit S als Pol (die in Formelsammlungen iibliche Darstellung): 

a 

p = 4a cosy - . 


How to trisect an angle with the trisectrix 

The angle § to be divided into three equal parts is 
drawn with its vertex in S and one ray on the x- 
axis. The other ray then intersects the trisectrix in 
P. 

ZSOP then is the third of the given angle (J>. 


cosy 
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Wahlt man S als Ursprung eines kartesischen Koordinatensystems (s. Zeichnung) so folgt mit 


p= -Jx 2 + y 2 


und cos \\j = — = 




x 2 + y 2 


rdie algebraische Gleichung fiir die 


Koordinaten des Kurvenpunktes P(x,y): 
4ax 


i- 


2 , 2 

x + y = 


V 


- a- 


2 , 2 

x + y 


x 2 2 

X — , also DERIVE: y 2 (a + x) - x 2 (3a - x) = 0 


Damit ist die Trisektrix von MacLaurin eine algebraische Kurve dritter Ordnung. 


Mit DERIVE lassen sich nun aus den Kurvengleichungen 

durch Auflosen der algebraischen Kurvengleichung nach y eine Parameterdarstellung der Kurve gewinnen 
(PLOT-Befehl!) oder eine andere Parameterform fmden 
die Asymptotik der Kurve untersuchen (x = -a) 

Nullstellen berechnen ((0/0) und (3 a/0)) 

mit Hilfe von S = — [r 2 dcp der von der Trisektrix eingeschlossene Flacheninhalt berechnen ( S = 3a 2 y[3) 

2 j 

Supported by DERIVE you may 

gain the parametric form of the curve (solve the equation for y) or find another one 
investigate the asymptote (x = a) 
find the zeros ((0/0) and (3a/0)) 
calculate the area enclosed by the curve. 


A DERIVE-Trisectrix-Session 

Expressions #1 through #8 are accompanied by a plot with slider bars for p and a. 


# 1 : 

* 2 : 

#3: 

/4 : 

*5: 


rayl := y = x-TAN(p) 

S := [2- a, 0] 

ray2 := y = (x - 2-a) -TANO-p) 

(SOLUTIONS (rayl a ray2, [x, y])) 

1 

2 

C0S(p) ■ (i - a - 4-a-SIN(p) ) a-SIN(p) ■ (2-C0S(2.p) + 1) 

. COS (p). COS (2- p) + SIN (p). SIN (2- p) ’ C0S(p) -C0S(2.p) + SIN(p) -SIN(2 ■ p) . 


* 6 : 


2 

C0S(t) ■ (i -a - 4.a.SIN(t) ) 

. COS (t). COS (2- t) + SIN(t) ■SIN(2-t) 


a-SIN(t) ■ (2-C0S(2-t) + 1) 

COS (t). COS (2- t) + SIN(t) -SIN(2 ■ t) . 


#7: Trigonoiietry := Collect 


a.COS(J.t) 

1 2 1 3 

a-SIN(3.t) 

- COS(t) 

C0S(t) 


This is another fine parameter form for the trisectrix. Now let's derive the analytic form by eliminating 
parameter t. 



D 
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#9 : Trigonometry := Expand 

#10: Trigpower := Cosines 


# 11 : 


# 12 : 


2 a -SIN (t) 

4-a.COSCt) - a, 4 -a -SIN(t) -COSCt) - 


cas(t) 


# 11 : 



Substitute for t in the second component of expression #1 1 : 

a -SIN(t) 

#14: y - 4.a-SIN(t)-CCS(t) - 


COS(t) 


#15: y = 4 -a -SIN 




ACOS 


■ COS 


ACOS 



and with little manipulation you can find the algebraic form. 


#16: y = 






#17: 


y = 


2 Y -o-a - x) 
#18: y = 


x + a 

The area 


# 20 : 


2 -n 



r a 



f * 

COS 

) 

J 

L h 

) } 


d4> 


= 3 ■ -/i ■ a 
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Kurzreferenz fur DERIVE User: 

Die Trisektrix von Maclaurin 



VECTOR (4 *a*COS (phi) -a/COS (phi) ,a, 0.5, 2. 5, 0.5) 


r=4*a*COS (phi) -a/COS (phi) 
[x, -x*SQRT ( 3a-x) / SQRT (x+a) ] 
y A 2 * (a+x) -x A 2 * (3*a-x) =0 


Polar form: 
Parameter form: 
Algebraic form: 


Ausblick: 

In der nachsten Folge des DERIVE-Kurven- 
Lexikons werden einige Konstruktionen vorgestellt, 
mit welchen die bisher vorgestellten Kurven 
einheitlich erzeugt werden konnen. Besondere 
Beachtung werden dabei einige Konstruktionen 
erfahren, die erst in den letzten Jahren unter Mithilfe 
der Computerprogramme DERIVE und Cabri- 
Geometre entdeckt wurden. 


Outlook: 

In the next part some constructions will be shown, 
which enable to produce the curves introduced till 
now in a homogenous way. We will give special 
attention to those constructions which have been 
found out in the last years with help of 
computerprograms as DERIVE and Cabri- 
Geometre. 


FlT^flT T Y FHt “Y _ 

|Rlg9bra|Calc|Oiher|PrgnIO|Cl0an Up| 


pfciean UpT~ 1 


1 y = x ■ tan(p) rayl y = tan(p)-x| 

1 y = (x - 2 ■ a) ■ tan(3 ■ p!) * ray2 

y = tan£3 ■ p) - (x - 2 ■ a)| 


■solyeCrayl and ray2, £x y}} 

~2 ■ a ■ cqs(p) ■ sin(3 ■ p) 

x - sin(p) ■ cos(.Z ■ p} - gos^p} ■ sin(3 ■ p) 

a ~2 ■ a ■ cqs(p) ■ sin(3 ■ p) ^ 

sinfpVcosf3-P^ - cos( pV sinf 3 ■ p1 


anc^ 


I MAIN FiliD AUTD FUHC 16/50 


u ~2 ■ a ■ sin(p) ■ sin(3 ■ p) ^ 

sin(p^ ■ gos( 3 ■ p} - gosCp} ■ sin(!3 ■ p} y 

~2 ■ a ■ s i nt» s i n( 3 ■ p!) 

sin£p> ■ cos(3 ■ p ) - gos(p) ■ sin£3 ■ p ) 


■ iCol lect(x') * x 

■ iCol leGt(y) + y 


MAIN RAD AUTD 


a ■ ( s i nC 4 ■ p^ + s i ntj 2 ■ p^) 
sin(2p) 

-a (gqsC4 ■ p!) ~ gqs{2 ■ p!)) 
sir\(2pf 


FUNC ifi/30 


FiTmoY fT^ "(~f^~"( fht y ft Y Y _ 'l 

|Rlgebra|CalG|Qihen|PngnIQ|ciean Up| 1 


1 y p * y 


a (sin(2 ■ t) + sinCtQ) 
sin(t) 

_ a (gos( 2 ■ t) - costt)) 
sin(!i ') 

1 tExpandtx!) 2 ■ a ■ gos(! t0 + a| 

1 solye(x_ = 2 ■ a ■ Gos(t) + a ? t) 

■-) and 


, (4 ■ En3 - l) n . |" a - x_ 

1= --sin- 1 !- 









D-N-L#13 


Thomas Weth: A Lexicon of Curves (3) 


pi 9 


On page 19 you can find the handheld derivation of the parameter form (which looks different as the 
others which have been found until now) and the implicit form. The next screen shots show parameter 
and polar form representations and the Tl-Cabri-produced trisectrix. 



J-LDTS H 


Plot Utei. n 
■■■'xt 1 =t. 

■/yU= tl 2 3} 

Jt + a 1 

■■■'xt2=t 

, yt2 =i^I| a = n 2 3> 

Jt + a 

xt3=B 


xt3<t >= 


MAIN RAD AUTD 

PAR 







9 

(}.S f 

I 

. . . . * 

28.61° 


MAIN -a DEG AUTD PDL 
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Another problem: (Dr. Weth) 


Given is the set M 


^ m n ^ 


y-n m j 


| ( m,n ) e C 2 


of 2 x 2 matrices. 


Wanted are all solutions of + 1=0 (with 1 = the 2 x 2 identity matrix). All my attempts to solve 
this equation with DERIVE failed. Also failed the attempt to solve the equivalent equation 
X + X’ 1 = 0. 


Calculating by hand I found solutions: 


^ icos a 

±sin a ^ 

und 

r i cos a 

±isin 

v + sin a 

-i cos a y 


±isin a 

-icos 


Maybe there is anybody among the DERIVE Users and - specialists all over the world who knows the 
solution. Many thanks for your help. 
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Computer Algebra and 
Spreadsheets (e.g. Excel) 

David Sjostrand, Onsala, Sweden 

In this paper we present some examples of how to use expressions, generated with a 
computer algebra system, DERIVE, in a spreadsheet program, MS Excel. These examples 
will appear in a coming book, Mathematics with Excel , (Studentlitteratur, Sweden). 


1. Systems of linear equations 

The system of linear equations 
a x x + b x y + c x z - d x 
< a 2 x + b 2 y + c 2 z - d 2 
a 3 x + b 3 y + c 3 z = d 3 

can easily be solved in DERIVE. 

The solutions will be expressions in a\, a 2 , d 3 , which is suitable if you take into account that Al, 

A2, ... and so on are cell references in a spreadsheet. 

Steps to be taken in DERIVE 

DERIVE should be in Word Input Mode. This can be achieved via Options > Mode Settings > Input. 
Edit the expression 

[al*x + bl*y+cl*z=dl,a2*x + b2*y+c2*z=d2 ,a3*x + b3*y+c3*z=d3] 

A system of linear equations is a vector of equations. Solve the system using the Solve-button 
(Looking Glass in the menu bar, highlight x, y and z .). The solution is given as a vector (list). Take 
care that the formulae for x, y and z will be separate expressions. You can do in the following way: 


#4: 


#5: 


# 6 : 



bl • (c2*d3 - 

- c3*d2) 

1 + b2* (c3 *dl 

- cl*d3) + 

b3 • (cl *d2 

: - c2*dl) 


al • (b2*c3 - 

- b3-c2j 

1 + a2*(b3*cl 

- bl*c3) + 

a3 • (bl *c2 

: - b2.cl) 

al 

•(c2*d3 - c: 

i *d2) + 

a2.(c3.dl - 

cl*d3) + a3 

• (cl *d2 - 

c2.dl) 

al 

♦(b2*c3 - b3*c2) + 

a2* (b3 *cl - 

bl*c3) + a3 

• (bl *c2 - 

b2*cl) 

al 

• (b2*d3 - b: 

• d2) + 

a2> (b3 *dl - 

bl*d3) + a3 

• (bl >d2 - 

b2*dl) 

al 

• (b2*c3 - b: 

■ c2) + 

a2* (b3 >cl - 

bl-c3) + a3 

• (bl >c2 - 

b2*clj jH 


bl • (c2*d3 - 

c3 *d2) 

+ b2*(c3*dl 

- cl *d3) + 

b3 • (cl *d2 

- c2.dl) 


al*(b2*c3 - 

b3 *c2) 

+ a2*(b3*cl 

- bl*c3) + 

a3 • (bl *c2 

- b2*cl) 


al • (c2*d3 

- c3*d2) + a2*(c3*dl - cl*d3) 

+ a3*(cl*d2 - c2*dl) 


al • (b2*c3 

- b3*c2 

0 + a2*(b3*c 

1 - bl *c3) 

+ a3*(bl*c 

2 - b2.cl) 


al • (b2*d3 - 

b3 *d2) 

+ a2.(b3.dl 

- bl*d3) + 

a3 * (bl *d2 

- b2*dl) 


# 7 : 


al * (b2*c3 - b3*c2) + a2*(b3*cl - bl*c3) + a3*(bl*c2 - b2*cl) 
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The expressions for x, y and z can be saved in BASIC-format and then be used in a spreadsheet. Save 
the formulae under the name LINEQU.BAS. (File > Write > Basic file) 




Save 


Expressions 

r State Variables 


C All 

W Expressions 


Selected 


v Make backup 


Netzwerkunngeb Dateiname: 
ung 

|linequ.bas| 

3 

Speichern 

Dateityp: 

| Basic file f.bas) 

"31 

Abbrechen 


m 


Steps to be taken in MS Excel 

(a) Open the file LINEQU.BAS. Cells Al, A3 and A5 should now contain the expressions for x, y 
and z as strings. 


Al - f* x=(b1 *(c2*d3-c3*d2)+b2*(c3*d1 - cl *d3)+b3*(d *d2- 


A 1 B C D E 

i 

x=(b1 *(c2*d3lc3*d2)+b2*(c3*d1 - cl *d3)+b3*(c1 *d2-c2*d1 Ma1*(b2*c3-b3*c2) 

2 


3 

Y=-(a1 *(c2*d3-c3*d2)+a2*(c3*d1 -cl *d3)+a3*fd *d2-c2*d1 ))/(a1 *fb2*c3-b3*c2 

4 


5 

z= fa 1 *fb2*d3- b3*d2) +a2*f b3*d 1 - b 1 *d3) +a3*fb 1 *d2- b2*d 1 ))/(a1 *fb2*c3-b3*c2) 

6 

□ 


(b) Move these expressions to the cells B7, B8 and B9 resepctively and remove the x, y and z but 
leave the = -sign to declare the expressions as formulae. 

(c) If we now, for example want to solve the system 
x + 16y + 3z = 6 

9x-9y + 2z = 4 
21x + 12y-35z = -2 

we have to enter the coefficients in the range Al :D3. 



A 

B 

C 

D 


1 

1 

16 

3 

6 


2 

9 

-9 

2 

4 


3 

21 

12 

-35 

-2 

- 

4 






(d) You should find the result in cells B7:B9 in decimal form. If you would like to have the result in 
exact form (as fractions) then reformat cells B7:B9 using the respective dialogue box: 


Zellen formatieren 


Zahlen 

Kategorie 


Ausrichliung || SchriFt || Rahmen || Muster || Schutz 


Standard 

Zahl 

Wahrung 

Buchhaltung 

Datum 

Uhrzeit 

Prozent 

Bruch 

Wissenschaft 

Text 

Sonderformat 


ieispiel 

2006/3447 


lyp: 


# ?????/????? 


# . ##0 ■ -* # . ##0 j _-* j a 

_- :+: #.##0 #.##0 J — ' 

_-* # . ## 0, 00 -* # . ## 0, 00 _-* 11 

#.## 0,00 #.## 0,00 
00000 

# ???/??? O 


BenutzerdeFiniert 
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(e) The : 

solution is: 


#8: S0LVE(x + 16-y + 3-z = 6 a 9-x - 9-y + 2-z = 4 

7 

x = 

2006/3447 

a 21-x + 12-y - 35-z = -2, [ x , y, z]) 

8 

y = 

850/3447 


9 

Z = 

188/383 

2006 850 188 

#9: x = a y = a z = 


3447 


3447 


383 


If you vary the entries in range A1 :D3 then you can solve every non singular system of linear 
equations consisting of three equations with 3 unknown. 


2. The circumcircle of a given triangle 

Let (xi, yi), (x 2 , y 2 ) and (x 3 , y 3 ) be three given points in the plane. Let (xM, yM) the center of the 

circumscribed circle of the triangle defined by the given points. 

xM andyM can be regarded as functions of the 6 variables x u y u x 2 , y 2 , x 3 andy 3 . 



Steps to be taken in DERIVE 

In DERIVE we can derive the expressions for xM andyM. 

As the three vertices of the triangle have equal distances to point (xM, yM) we see that we can get the 
desired position if we solve the following system of equations, which is in fact linear in xM andyM. 


2 2 2 2 

#10: equl := (xl - xin) + (yl - yin) = (x2 - xin) + (y2 - yin) 

2 2 2 2 

#11: equ2 := (xl - xin) + (yl - yin) = (x3 - xin) + (y3 - yin) 

(SOLUTIONS (equl a equ2, [xin, yii])) 

# 12 : 1 


Use SOLVE or SOLUTIONS to solve the system with respect to xM andyM You will receive two 
bulky expressions for xM andyM. If you apply the EXPAND command on the numerator and 
denominator, they will look like this: 


#18: xii 


22222222 222 2 
xl -y2 - xl -y3 - x2 -yl + x2 -y3 + x3 *yl - x3 -y2 + yl -y2 - yl -y3 - yl-y2 + yl*y3 + y2 *y3 - y2-y3 


2 • xl ■ y2 - 2xly3 - 2-x2.yl + 2x2y3 + 2-x3.yl - 2-x3-y2 

22 22222 22222 
xl -x2 - xl -x3 - xl-x2 + xl-x3 - xl-y2 + xl«y3 + x2 *x3 - x2-x3 + x2-yl - x2-y3 - x3-yl + x3*y2 


#19: yii = — 


2 ■ xl • y2 - 2-xl.y3 - 2-x2yl + 2-x2y3 + 2-x3.yl - 2-x3.y2 
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Select (Highlight) expressions #18 and #19 and save them in BASIC format under the name, say 
CIRCUMS.BAS. 


Steps to be taken in MS Excel 

Open the file CIRCUMS.BAS. Cells A1 and A3 are containing the expressions for xM and yM as 
strings. These strings can again easily be converted to an MS Excel formula simply by removing the 
“xm” and “ym”. 

: U j J | ■ | ^ - | ■ Bearbeitung zurucksende 

A3 - £■ =-(xi^*x2-xi^2*x3-xra^+xra^-xry; 

Then you will receive an error message, 
because you have nothing but empty cells in 
the range XI : Y3. Now you can . . . 



A I 

B Y3+2*X3*Y 1 -2*X3*Y2) 

1 

' #DIV/0! 





2 






3 

' #DIV/0! 1 





4 

□ 






Create a worksheet to plot a triangle with given vertices together with its circumcircle. 

1. Move the expression for yM to cell B3 and the expression for xM to cell A3. You can drag with 
the mouse or use the commands Edit Cut and Edit Paste. Enter xM and yM in the cells A2 and 
B2 respectively. 

2. Plot a triangle with vertices (3,1), (9,4) and (2,6). Enter the coordinates in the range XI : Y3 and 
the formula =X1 and =Y1 in the cells X4 and Y4 (to close the triangle, similar as in DERIVE). 
Then move the contents of range XI :Y4 to the range C2 : D5 by again applying Cut and Paste. 

Select the range C2 : D5 and plot the triangle as a (XY)-diagram. 


3. We plot the circle as another (XY)-diagram and have to provide the points forming the 

circumference. The points will be defined by the parameter representation of the circle. The 
equation of the circle is: 


(x = xM + R cos t 
\y = yM + Rsint 


with R = yxl-xMf +(yl-yM f 


Enter R in cell 7 and the formula =SQRT((C2-A3) A 2+(D2-B3) A 2) in cell C8. (The German 
version of MS Excel needs WURZEL(...) instead of SQRT.) 

Enter Circumscribed Circle in cell A10 and t, x and y in the cells A1 1, B1 1 and Cl 1 
respectively. Then enter 0 in cell A12 and the formula = PI()/50+A12 in cell A13 and copy it to 
the range A13 : A1 12 (= 2n). 

Then enter the formulae =$A$3+$C$8*COS(A12) and =$B$3+$C$8*SIN(A12) in cells B12 
and C12 respectively. Copy the contents of B 12 : C 12 to the range B1 12 : Cl 12 and select range 
B12 : Cl 12. Move the cursor - with the Ctrl-key pressed - to the top right corner of the selected 
range. The cursor changes to a little Now you can drag the table into your diagram and add 
the points of the circle. If you change any of the coordinates then the triangle and subsequently 
the circle will change accordingly. You might find a superfluous segment connecting the first 
point of the circle with the last point of the triangle. Then select this segment and make it 
invisible by reformatting. 


Now you should have something like this: 
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Adjust scales manually in order to receive a figure where the circle really looks like a circle. 
Select the axes and give the command Format Scale. 


B12 - ft = $ A$3 +$ C $8*C 0 S (A1 2) 



A 

B 

C 

D 

E 

F 

G 

H 


1 



x-Coord 

y-Coord 






2 

xM 

yM 

3 

1 






3 

5,29411765 

4,05882353 

9 

5 






4 



2 

6 






5 



3 

1 






6 










7 



R 







8 



3,82352941 







9 



1 







10 

Circumscribed Circle 



- 

11 

t 

x v 


12 

o| 

1 9.117647061 4.05882353 

9 n 


13 

0,06283185 

9,1101022 4,29890493 


14 

0,12566371 

9,08749739 4,53803883 

8 - 

15 

0,18849556 

9,04992184 4,7752815 


16 

0,25132741 

8,99752385 5,00969663 

7 - 

17 

0,31415927 

8,93051021 5,2403591 

6 - 

18 

0,37699112 

8,84914539 5,46635858 


19 

0,43982297 

8,75375049 5,68680317 

5 - 

20 

0,50265482 

8,64470201 5,90082287 


21 

0,56548668 

8,5224303 6,10757304 

4 - 

22 

0,62831853 

8,38741792 6,30623773 


23 

0,69115038 

8,24019769 6,4960329 

3 - 

24 

0,75398224 

8,08135063 6,67620952 

2 - 

25 

0,81681409 

7,91150364 6,84605652 


26 

0,87964594 

7,73132702 7,00490358 

1 - 

27 

0,9424778 

7,54153185 7,1521238 





28 

1 ,00530965 

7,34286716 7,28713619 

n 






29 

1,0681415 

7,13611699 7,40940789 

C 

)1 23456789 10 

30 

1,13097336 

6,92209729 7,51845638 



We can now reap the fruits of our efforts, You can plot a new triangle and its circumscribed circle 
just by entering new coordinates in the cells in the range C2 : D4. 



The DERIVE slider bar version from 2007 
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3. The incircle of a given triangle 

What we did in 2. with a given triangle and its circumscribed circle, you can do in a corresponding 
way with a given triangle and its inscribed circle. 

Let (xi, yi), (x 2 , yi) and (x 3 , y 3 ) be three given points in the plane. Let (xM, yM) the center of the 

inscribed circle of the triangle defined by the given points. 

xM and yM can be regarded as functions of the 6 variables x\,yu * 2 , yi, x 3 andy 3 . 


Steps to be taken in DERIVE 

In DERIVE we can derive the expressions for xM and yM. 



The side opposite to point (xi,yi) has length /i, i = 1, 2, and 3. The intersection point between the 
bisector from (xl,yl) and the opposite side is (x4,y4) and the intersection between the bisector from 
(x2,y2) and the opposite side is (x5,y5). The centre of the incircle is the intersection point of the 
bisectors of the triangle. 

According to the theorm of bisectors point (x4,y4) divides side 11 connecting (xl,y2) with (x2,y2) in 
the ration 13/12. Therefor we get the equations 

#25: [12- (x2 - x4) = 1i-(x4 - xi) , 1 2 ■ (y2 - y4) = 1i-(y4 - yi)] 


They form a system of linear equations in x4 and y4. It can easily be solved: 

*26: S0LVE([12 ■ (x2 - x4) = 1i-(x4 - xi) , 12- (y2 - y4) = 1i-(y4 - yi)], [x4, y4]) 


#27: 


!2-x2 + 1 i ■ xi !2-y2 + 1i-yi 

x4 = a y4 = 


12 + 1i 12 + 1i 

In a corresponding way we get a system of equations for x5 andy5. If you solve it you will get: 

*i0: 


11-xl + 1 i ■ xi 11-yl + 1i-yi 

x5 = a y5 = 


11 + 1i 


11 + 1i 


Enter the following assignments for the coordinates: 

!2-y2 + li-yi 


*il: 


!2-x2 + 1 i ■ xi 


x4 := 


y4 := 


11-xl + li-xi 11-yl + li-yi 

xS := , y5 := 


12 + 1i 12 + li 11 + li 11 + li 

and set up the next system of equations formed by the equations of the two bisectors: 
y4 - yl y5 - y2 


f 33 : 


y - yl 


x4 - xl 


Cst - Xl), y - y2 


x5 - x2 


(x - x2) 
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#34: SOLVE 


#35: 


(b 


yi 


y4 - yl 


(x - xl), 


x4 - xl 

11-xl + 12 ■ x2 + 13-x3 


11 + 12 + 13 

Rename the solutions as xM and yM 


* y 


ys - y2 ] 

y - y2 = (x - x2) , 

x5 - x2 

1 1 ■ yl + 1 2 ■ y2 + 13-y3 ' 


11 + 12 + 13 


[x, y] 


11- xl + 12- x2 + 13- x3 

#36: xiii = 

11 + 12 + 13 

11- yl + 12- y2 + 13-y3 

#37: ym = 

11 + 12 + 13 

and finally substitute for the lengths of the sides /l, 12 and 13. 

/1 = V(x3-x2) 2 +(j3-j2) 2 , 12 = -y/(jc3 - xl) 2 + (j3 - jl) 2 , 13 = y[{x2 - xl) 2 + (yl - jl) 2 


#38: xm 


2 2 2 2 2 2 
V((x3 - x2) + (y3 - y2) ).xl + 7((x3 - xl) + Cy3 - yl) ).x2 + V((xl - x2) + (yl - y2) )-x3 


2 2 2 2 2 2 
7((x3 - x2) + (y3 - y2) ) + V((x3 - xl) + (y3 - yl) ) + V((xl - x2) + (yl - y2) ) 


#39: 


ym r 


2 2 2 2 2 2 
7((x3 - x2) + (y3 - y2) ).yl + ^((x3 - xl) + (y3 - yl) ).y2 + V((xl - x2) + (yl - y2) ).y3 


2 2 2 2 2 2 
V((x3 - x2) + (y3 - y2) ) + V((x3 - xl) + (y3 - yl) ) + V((xl - x2) + (yl - y2) ) 


Save the expressions in BASIC format as INCIRCLE.BAS. 

Steps to be taken in MS Excel 

Open the file INCIRCLE.BAS and do as before. SQR must be replaced by SQRT or WURZEL. 


For the radius of the incircle we 
can use the formula 

A 

r = — 

P 


where 

^ = V J p(p-/l)(^-/2)(p-/3) 
and 


/1 + /2 + /3 

p- 

2 

is half the perimeter of the 
triangle. 


A3 - J* = (WU RZE L((C4- C3) A 2 +(D4- D3) A 2)*C2 +WU RZE L((C4- C2) A 2 +(D4- D2) A ™*C3 +WU F 


A | B WU RZE L((C4- C2) A 2 +(D4- D2) A 2) +WU RZE L((C2- C3) A 2 +(D2- D3) A 2)) 


1 



x-Coord 

y-Coord 



2 

xM 

yM 

1 

-1 





3 

3.7499954 

2,4218906 

9 

2 





4 



2 

6 





5 



1 

-1 





6 

11 = 

8,0622577 







7 

12 = 

7,0710678 

R = 

2,23842813 





8 

13 = 

8,5440037 







9 




7 _ 





10 

P = 

11,838665 







11 




6 - 





12 

INCIRCLE 



5 - 





13 

t 

X 

y 






14 

0 

5,9884235 

2,4218906 

4 - 





15 

0,0785398 

5,9815232 

2,5975156 

3 - 





16 

0,1570796 

5,9608647 

2,7720579 






17 

0,2356194 

5,9265756 

2,9444412 

2 - 





18 

0,3141593 

5,878867 

3,1136029 

1 - 





19 

0,3926991 

5,8180333 

3,2784999 






20 

0,4712389 

5,7444495 

3,4381157 

0 - 



7 8 9 


] 

21 

0,5497787 

5,6585691 

3,591466 

-1 - 

V 

-J H 

► O O 

1 

22 

0,6283185 

5,5609218 

3,7376056 






23 

0,7068583 

5,4521095 

3,8756333 

-2 - 
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4. Systems of two non linear equations with two unknown 

A system of equations canbe solved numerically with the method of Newton Raphson. For the system 

I7(*.j0 = o 
tg(*>y) = o 

consisting of two equations with two unknown, x and y, and the solution ( a. /:!) w e get according to 
Newton Raphson's method that (x n , v n ) — > (a, /?) for n — » co if the point (xi, Vi ) is sufficiently close to 
(a,/3). 

x n and y n , n > 1 are given by the formulae 

f(x„,y„)gy'(x„,yJ-g(x n ,y B )f y '(x B ,y n ) 

f x '(x„,y n )g y '(x n ,y„)-g x '(x„,y n )f y '(x n ,y„) 

= f x '(x n ,y n )g(x n ,y n )-g x '(x n ,y n )f(x n ,y n ) 
f x , (x n ,y n )g y '(x n ,y n )-g x , (x n ,yj f y '(x n ,y n ) 


Steps to be taken in DERIVE 

Enter the expressions 

d d 

f g - g f 

d yl d yl 

xl - 

( d ^ d ( d ^ d 

f g - g f 

Ld xl J d yl Id xl J d yl 



You can enter the expression f simply by typing DI F (f , xl) . 

Id xl J 

Suppose that we want to find approximate solutions to the system (User Manual DERIVE, page 143) 

J 3xy-2x 2 +4siny + 6 = 0 
[3x 2 -2 xy 2 +3cosx + 4 = 0 

Enter the expressions 

[ 2 2 2 1 
Lf := 3-xl-yl - 2-xl + 4-SIN(yl) + 6, g := 3-xl - 2-xl-yl + 3-C0S(xl) + 4j 

and simplify the expressions from 

above. You will get formulas for x 2 and 3.cos(xi).(4-tosCyi) + s-xi) + i.xi.siN(xi).(4.cosCyi) + 

0 • *45: 

y 2 . Save these two expressions as 

NEWRAPH.BAS and start MS Excel. 3 -sin(*d ■ C4.costyi) + 3-xi) ■ 


Part of the simplified value for x 2 : 


■yl - 9) + 12 ■ (2 -yl + 1)) 
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Steps to be taken in MS Excel 

Open the file NEWRAPH.BAS. Letxi = 3 and y\ = 2 in cells X2 and Y2. Enter the expressions forx 2 
andj 2 ’ in cells X2 and Y2. (Don't forge to add the leading =-sign.) 

Calculate the sequence of x- andy-values by copying the contents of range X2 : Y2 to the range 
X2 : Y10. 


We can see that our initial values xi = 3 and y \ = 2 yield an approximate solution: 


x « 4.9543, y « 2.8120 

If you vary the content of the cells in XI and Y 1 you can find other solutions if there are any. 


I copied the cells into the V and W-columns 
and found a second solution. 

Then I plotted both equations together with 
the found solution points. 


m 

)+4*(4-3*V1 A 2)*COS(W1)+V1 *(1 6*W1 *SIN(W1 )+V1 A 2*(8*W1 -9 
l*(8*W1-9)-3*W1 A 2))) 


1 

-1 

3 

2 

I -0,2251 G987f-2 ,20326553 

6,53694238 

3,57925193 

0,75731861 

-2,67973572 

5,28841854 

3,06177454 

0,48793859 

r -2 ,707331 36 

4,98122898 

2,83980263 

0,48563645 

r -2,7530482 

4,95455164 

2,81238333 

0,48573914 

r -2 ,752621 39 

4,95424665 

2,81204193 

0,48573916 

r -2 ,752621 39 

4,9542466 

2,81204188 

0,48573916 

r -2 ,752621 39 

4,9542466 

2,81204188 







Just for fun I solved the system using the 
Excel Solver and then using the tools 
provided by DERIVE. 

Josef 

NEtfTONS([f, g], [xl, yl], [3, 2], 5) 

3 2 

6 . 536942383 3.579251930 

5.288418538 3.061774538 

4 . 981228980 2 . 839802628 

4.954551640 2.812383328 


: _J 'Id ■- - (|_!l S Bearbeitunc 


A2 - f* =3*A5*B5-2*A5 A 2 +4*5 1 N (B5) 



A | 

B 

C 

D 

1 





2 

-6,000000021 

1 -6 



3 

-4 ,0000001 ' 

-4 



4 

X 

V 



5 

4,95424661 

2,81204188 


1 ,1 042E-1 4 


FIRST (REVERSE(NEWT0NS([f , g], [xl, yl], [i, -2]))) 
[0.4S57B9161S, -2 . 7526213B6] 


4 . 954 24664 S 2 . S12041929 
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Minimization of a "Flat Function" 

A Classroom Experiment with DERIVE 

Jose M. Cardia Lopes, Porto, Portugal 

Edgar and Himmelblau 1) present a formulation problem where the operation cost, C, of 
a hypothetical chemical plant is a function of the amount of product per batch, P. When the 
problem is solved we obtain the following function: 

c _ 400, 000 + 8QOP 07 + 75,000, 000 ( '14 + 2_P ft4 J 


With this function we can motivate our students for th need and particularities of some 
common numerical approaches. And DERIVE is a good help to do it in the classroom. 

Suppose we want to claculate the value of P that minimizes the operating cost, C. The first 
idea is to determine the derivative, dC/dP, and solve it to zero: 

0.4 

0.7 15000000- (14 + 2-p ) 

#1: c := 400000 + S00-p + 


In 1994 C. Lopes asked for “zeroing" the derivative in MIXED & EXACT Mode. 


1 2: — c 

dp 




17/10 2/5 

SO- (7-p - 225000- p - 2625000) 


#4: Precision := Exact 

#5: Notation := Rational 


#6 : SOLVE 


17/10 2/5 

SO - (7 - p - 225000- p - 2625000) 


17/10 


2/5 


*7: 

7 ■ P 

- 225000- p 

= 2625000 v p = unit_circle ■« 

#8: 

Precision nixed 




#9: 

Notation := Flixed 






' 17/10 

2/5 



#10: 

SOLVE 

80- (7-p 

225000- p - 2625000) 




2 

* P 







P 



#11: 

p = unit_circle ■« ¥ p 

= 3856.931625 





17/10 

2/5 



#12: 

SOLVE 

SO- (7-p 

225000- p - 2625000) 


Real 


2 

i P> 






-V 

P 


j 


17/10 2/5 




#13: 

7-p 

- 225000- p 

= 2625000 v p = +» 




D 


Edgar & Himmelblau, Optimization of a Chemical Process, McGraw-Hill, Inti Edition, p28, problem 1.7 
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Surprise: DERIVE is unable to 
solve this equation (at the Exact or at the 
Mixed Mode)! But the function really has 
a minimum. In fact if we graph this 
function (and the first problem for most 
students, is to choose an appropriate 
scale) we can see that the minimum 
exists and it is close to P = 3850. Then: 
Why is it “impossible” to determine its 
mini-mum? To answer this question it is a 
good idea to plot dC/dP. 



S 10 A 5 


■6 10 A 5 
•4 10 A 5 
■2 10 A 5 

Derivative dC/dP 

T 1000 2000 5000 4000 5000 

-2 10 a 5 


As we can see in the figure the graph of dC/dP is nearly coincident with the axis (that is why 
the intersection is so difficult to find). And it is the right time to speak a little about ill 
conditioned problems ... 


Now we can try to use an interval 
method to solve this equation - 
we work in approximate mode (or 
we use NSOLVE). 


#14: 


NSOLVE 


17/10 

80-C7-P 


2/5 

225000- p - 2625000) 

2 

P 


p, Real 




t 15: p = 3856. 931625 


Now the students are (probaby) motivated to discuss the particularities of the two classes of 
numerical methods for zeroing a function: the iterative and the intervalar methods. 

We can also try another approach: to approximate the function by a Taylor polynomial (which 
is the “right” degree? and the “right” base point?), determine the derivative and find its 


zero(s). 


#20: 

TAYLORCc, p, 3850, 2) 



2 

6 

#21: 

0. 008704314509 -p - 67 

. 14356204 -p + 1.054558048-10 


d 2 

6 

#22: 

— (0. 008704314509- p 
dp 

- 67 . 14356204 - p + 1.054558048-10 ) 

#23: 

SOLVE (0 . 01740862901 ■ p 

- 67.14356203, p, Real) = ll.MLm.llfl HUM) 

#24: 

TAYLORCc, p, 3850, 3) 

r d 


#25: 

SOLVE — TAYLORCc, p, 
Idp 

3850, 3), p 

#26: 

p = 3856.931669 v p = 

6373.566236 



Perhaps the Taylor approximation is not the best approach to this problem but it is always a 
possible approach and a new and enriching topic to discussion in the classroom. 



D-N-L#13 


Johann Wiesenbauer: Titbits (1) 


p31 



Johann Wiesenbauer, Vienna, Austria 


This is the first part of a lecture held by Dr. JohannWiesenbauer from the Technical University, 
Vienna, in the frame of the 3rd International Conference on School Mathematics in Vienna, 
23.2. - 25.2.1994. Dr. Wiesenbauer showed some examples of sophisticated ’’programming” with 
DERIVE. The DERIVE functions should be loaded as Utility files. In the next issue you will find the 
second part of this lecture accompanied by Dr. Wiesenbauer 's explanation and interpretation. You will 
marvel how he is handling and investigating - supported by DERIVE - huge numbers with ease. The 
word has Dr. Wiesenbauer: 

Titbits(l) - then and now 

Johann Wiesenbauer, Vienna University of technology, February 2007 

This is a revised version of my Titbits(l) in the DNL #13. Note that if there are two functions 
with the same name, the f irst one is usually the old version (sometimes with small notational 
changes), which is then overwritten by an updated version that makes full use of all the features 
of Df W 6.10 along with other improvements. 

dellast(v) returns a copy of the vector v without its last component 

deDast(v) := VECTOR (v , k, 1, DIH(v) - 1) 

#1: k 

*2: deDast(v) := DELETE (v, -1) 

euclid(a,b) returns the list of all equations in the Euclidean algorithm. 

eucl-id(a, b) := dellast (ITERATES (IF(x = 0, x, lx , x , FLOOR(x , x ), x , x , HO0(x , x )]), x, [a, =, FLOOR (a, b) , 
#3: 6 L 4 2 4 6 6 5 4 6 J 

b, +, ROD (a , b)]» 

#4: euc lid (1234567891 , 234567891) 

123-4567891 = 5 254567891 + 61728456 

254567891 = 5 61728456 + 49582585 

61728456 = 1 49582585 + 12545855 

49582585 = 5 12545855 + 12545024 

12545855 = 1 12545024 + 829 

12545024 = 14891 829 + 585 

*5: 

829 = 2 585 + 59 

585 = 6 59 + 51 

59 = 1 51+28 

51 = 1 28 + 5 

28 = 9 5 + 1 

5=5 1 + 0 
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euclid(a, b, q_, r_, s_ := [], t_) := 

Loop 

If b = 0 

RETURN REVERSE (s_) 
q_ := FLOORCa, b) 
r_ := HOD (a , b) 

*6: If r_ = 0 

t_ := [a, q_, b, MM , ""] 

t_ := [a, q_, b, r_] 

s_ := AD JOIN (t_, s_) 
a := b 
b := r_ 

*7: euclid(1234567S91, 234567891) 


1234567891 = 

5 

■ 234567891 

+ 

61728436 

234567891 = 

3 

■ 61728436 

+ 

49382583 

61728436 = 

1 

■ 49382583 

+ 

12345853 

49382583 = 

3 

■ 12345853 

+ 

12345024 

12345853 = 

1 

■ 12345024 

+ 

829 

12345024 = 

14891 

829 

+ 

385 

829 

2 

385 

+ 

59 


fibpair(s) returns the largest pair of adjacent Fibonacci numbers not greater than s. (Here the 
update is rather an alternative than an improvement!) 

f-ibpa-ir(s) := ITERATE ( I F(ir + y > s, y, [y + y , y T ) , y, [1 , 0]) 

* 9 : 12 L 1 2 lj 

100 

#10: fibpa-irClO ) 

#11: [9216845717656874712980450562726202415567360565980794777111390850331644813674856981646960226192287360, 

5696323922575865414847061494575945648081290145228607189038829076215134884313127297923138542545712321] 

f-ibpanrCs, f_ := [1, 0]) := 

Loop 

#12: IF I(F_) > s 

RETURN f_ 

f_ := [1, i; 1, 0]-f_ 

100 

#13: fnb|>a-irClO ) 

#14: [9216845717656874712980450562726202415567360565980794777111390850331644813674856981646960226192287360, 

5696323922575865414847061494575945648081290145228607189038829076215134884313127297923138542545712321] 

In spite of the now built-in function fibonacci(n), the following two (original) ways of computing 
Fibonacci numbers are still interesting for didactic reasons, as they show how the "golden 
section" number X = (l+75)/2 comes into play. 

f-ibl(n) := 

75 

f 1 f 1 + 75 Tn 

*16: fib2(n) := FLOOR +0.5 

L 75 L 2 J 

*17: 

*18: 


* 19 : 


f-ibl(lOO) = 354224848179261915075 
fib2(100) = 354224848179261915075 
FIBONACCI (100) = 354224848179261915075 
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eea (a , b) := del las t (ITERATES (IF(ROO (x , x ) = 0, x, [x , N0D(x , x ), x , x - FL00R(x , x )>x , x , x - FLOOR(x , 


1 2 4 6 5 


x )x 1), x, [a, b, 1, 0, 0, 1])) 

2 6j 

#30: eea(1234567 , 234567) 


1234567 

234567 

1 

0 

0 

1 

234567 

61732 

0 

1 

1 

-5 

61732 

49371 

1 

-3 

-5 

16 

49371 

12361 

-3 

4 

16 

-21 

12361 

12288 

4 

-15 

-21 

79 

12288 

73 

-15 

19 

79 

-100 

73 

24 

19 

-3207 

-100 

16879 

24 

1 

-3207 

9640 

16879 

-50737 


Well, what should I say more? It still works, which shouldn't be taken for granted in view of all 
the versions of Derive we have seen since version 2.56 used then! Needless to say that this 
piece of code would look completely different now, maybe like this. 

eea(a, b, q_, r_, s_ := []) := 

Prog 

a := [a, [1, 0]] 
b := [b, [0, 1]] 

Loop 

If FIRST(b) = 0 

*32: RETURN REVERSE (s_) 

s_ := AD JOIN ([FIRST (a), FIRST(b) , FIRST (REST (a)) , FIRST(REST(b)) ] , s_) 
q_ := FLOOR (FIRST (a) , FIRST(b)) 
r_ := a - q_-b 


*33: eea (1234567 , 234567) 


1234567 

234567 

[1, 0] 

[0, 1] 

234567 

61732 

[0, 1] 

[1, -5] 

61732 

49371 

[1. -5] 

[-3, 16] 

49371 

12361 

[-5, 16] 

[4, -21] 

12361 

12288 

[4, -21] 

[-15, 79] 

12288 

73 

[-15, 79] 

[19, -100] 

73 

24 

[19, -100] 

[-3207, 16879] 

24 

1 

[-3207, 16879] 

[9640, -50737] 


You will certainly have noticed that we introduced a new order of the colums as well as some 
pairings of numbers. This has the following meaning: At any time numbers in the first and second 
columns are the inner product [x,y][a,b] where [x,y] can be seen in the third and fourth columns, 
respectively. In particular, the gcd(a,b), which is 1 in this example, can be written as 

*35: [9640, -50737] -[1234567, 234567] = 1 

By the way, if you are only interested in this vector [9640, -50737], which is quite often the 
case in the applications, then you could also use the built-in function extended_gcd(a,b). It 
should also work for Gaussian integers, but sometimes it fails. This is a different story though, 
which I will tell you another time! 

*36: EXTENDED.GCD (1234567, 234567) = [1, [9640, -50737]] 
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Many, many years ago, I found an interesting article in the Scientific American about Public Key 
Cryptology. Specially the RSA- and the Merkle-Hellman- Algorithm were described. I was fascinated 
of this number theory application and wrote a "GWBASIC" program. (That was 1978 or 79 and I had 
time enough for programming, because there was no DERIVE and no DUG! ! !). For Merkle-Hellman 's 
method I needed the EEA-algorithm applied on two relative prime numbers. Having attended 
Dr. Wiesenbauer's lecture I tried to find my program "KRYPTO.BAS". Here you can see some lines 
BASIC-Code (the Berlekamp!). The full program is very nice and demonstrates in an easy to 
understand way the method how to produce a public key, how the system works, how to encrypt and 
to decrypt a message. If anybody is interested in this program, then please let me know. Josef 

REM A bit adapted for QBASIC 
DEFLNG A-Z 

INPUT "Two numbers with GCD =1: ", Zl, Z2 

XO = 1: Y1 = 1 : 

A = Zl: B = Z2 
start : 

Q = INT (A / B) : R = INT (A - Q * B + .5) 

X2 = XO - Q * XI: Y2 = YO - Q * Y1 

IF Y1 > 0 THEN 01$ = "" ELSE 01$ = "-" 

IF Y2 > 0 THEN 02$ = "+" ELSE 02$ = "-" 

IF R < 2 THEN PRINT R; " = " ; 01 $ ; ABS (X2 ) ; "x" ; Z 1 ; 02 $ ; ABS ( Y2 ) ; "x" ; Z2 : END 
A = B: B = R: XO = XI : YO = Y1 : XI = X2 : Y1 = Y2 : GOTO start 


Two numbers with GCD = 1: 1234567,234567 
1 = 1234567 x 9640 - 234567 x 50737 


The BASIC-program was my contribution in 
1994. In these days when the CAS-calculators 
are very welcome tools in math education 
(unfortunately not in all schools and all countries) 
I can offer the Extended Euklid Algorithm on the 
CAS devices TI-92 / Voyage 200. Josef 


|Fi7»ir|Y F£T Y FjT Y FHt Y FE Y Ffi T T "l 

| t |fl 1 gebra |Ca 1 c lother |Pr gm 1 0 |c 1 ean Up | 

■ eeaCl 234567, 234567) 

Cl 9640 1234567 "50737 

234567 > 

■ eeatl 234563, 234563) 

£3 "19 123456S 100 

234563 > 

■ gcd( 1234563 ? 234568) 

3 

■ - 19 ■ 1234568 + 100 ■ 234568 

8 

"19*1234568+100*234568 


CRYPTO RAD EKACT FUNC H/30 



Message 2997 From SOFT WAREHOUSE to JERRY GLYNN about GARBAGE COLLECTIONS 

You asked why DERIVE recycles memory (i.e. garbage collect) even though the percentage of free memory 
displayed on the bottom of the DERIVE screen never goes below 99% for a particular problem? 

When DERIVE is simplifying a mathematical expression, it generates a lot of data stmctures that are used only 
temporarily. Space is allocated for these stmctures as they are generated until the end of free memory is 
reached. At that point a garbage collection occurs to collect these no longer needed data stmctures. The two- 
pass, compacting garbage collector begins by marking all stmctures that are still being used by the program (i.e. 
those pointed to by program variables or the stack). Next the garbage collector makes a linear sweep through 
memory collecting all stmctures that are not marked (i.e. garbage) and makes the space available for reuse. 

At this point, * after* the garbage collection, DERIVE computes and displays the percentage of free memory 
statistic. Note that this is the *only* time a meaningful statistic can be calculated since only then is the size of 
the active data stmctures known. 
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Therefore if a lot of data structures are being produced, many garbage collections will occur. But, if most or all 
of the data structures being produced are temporary, the garbage collection will free up most of memory and the 
percent free will remain high. And, it should remain high since the system is *not* in danger of running out of 
memory. 

Aloha, A1 Rich, Soft Warehouse, Inc. 

DNL: Answers to inquiries of Mr Herdt and Mr Lopes: 

In one of the later messages in the BBS I found the similar problem: 

Message 3374 From SOFT WAREHOUSE to KEITH WILLIAMS about TRANPOSE OPERATOR BUG 

Thank you for your fax letter dated December 24, 1993 concerning a problem with DERIVE's transpose 
operator. The problem arose because DERIVE version 2.59 improperly transposes ragged arrays. Ragged arrays 
are vectors of vectors not all having the same number of elements. The problem did not arise in earlier versions 
of DERIVE because those versions did not even attempt to transpose ragged arrays. 

The problem has been resolved and will be included in the next version of DERIVE. Please send me your 
mailing address and I will send you a complementary update as soon as it is available. 

Aloha, A1 Rich, Soft Warehouse, Inc. 

I sent a FAX to SWHH in the evening and thanks modem times same night at 0322 am I received an answer 
from A1 Rich: 

In response to requests by users, beginning with version 2.59 DERIVE transposes a vector to an n by 1 column 
matrix. The transpose of an n by 1 column matrix is a n by 1 row matrix. Since a 1 by n row matrix is not 
equivalent to an n element vector they cannot be subtracted and their difference is not the zero vector. 

Perhaps DERIVE should transform the vector into a row matrix so the subtraction can occur. However, this is 
not valid if the vector has symbolic elements that couls themselves be replaced by vectors. 

The screen problem that occurs when mnning DERIVE under Soft-PC must be a bug in Soft-PC, not DERIVE. 
Apparently Soft-PC does not correctly emulate the screen services provided by the IBM-PC BIOS. I recommend 
Mr. Lopes contact the authors of Soft-PC to see if they can resolve the problem. Sicerely, A1 Rich 


Albert Floch'lay, Dirinon, France 


Je profite de la presente lettre pour vous remercier d'avoir envoy ez ma lettre du 6 Juillet 1993 au Dr Kutzler. 

Celui-ci a tres bien compris Lobjet de ma mauvaise humeur, mai la prix excessif de la mise a jour n'etant pas du 
a SOFT WAREHOUSE mais a Limporteur francais, je continue a faire confiance a DERIVE (tout en changeant 
fournisseur). Avec mes meilleurs voeux pour lannee 1994. 



The DERIVE Owl of Wisdom 






